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ABSTRACT 

Absolute ages of young stars are important for many issues in pre-main sequence stellar and cir- 
cumstellar evolution but are long recognized as difficult to derive and calibrate. In this paper, we 
use literature spectral types and photometry to construct empirical isochrones in HR diagrams for 
low-mass stars and brown dwarfs in the rj Cha, e Cha, and TW Hya Associations and the /3 Pic 
and Tuc-Hor Moving Groups. A successful theory of pre-main sequence evolution should match the 
shapes of the stellar loci for these groups of young stars. However, when comparing the combined 
empirical isochrones to isochrones predicted from evolutionary models, discrepancies lead to a spectral 
type (mass) dependence in stellar age estimates. Improved prescriptions for convection and boundary 
conditions in the latest models of pre-main sequence models lead to a significantly improved correspon¬ 
dence between empirical and model isochrones, with small offsets at low temperatures that may be 
explained by observational uncertainties or by model limitations. Independent of model predictions, 
linear fits to combined stellar loci of these regions provide a simple empirical method to order clusters 
by luminosity with a reduced dependence on spectral type. Age estimates calculated from various sets 
of modern models that reproduce Li depletion boundary ages of the /3 Pic Moving Group also imply 
a ~ 4 Myr age for the low mass members of the Upper Sco OB Association, which is younger than 
the 11 Myr age that has been recently estimated for intermediate mass members. 

Subject headings: stars: pre-main sequence — stars: planetary systems: protoplanetary disks matter 
— stars: low-mass 


1. INTRODUCTION 

Age estimates for pre-main sequence stars are critical 
in assessing the timescales for the formation of stellar sys¬ 
tems and the evolution of molecular clouds. The fraction 
of embedded sources determines how long the dominant 
phase of stellar growth lasts and how quickly the envelope 
disperses. Similarly, when combined with disk fractions, 
age measurements place strict constraints on the time 
available for giant planets to form. Ages are also impor¬ 
tant in quantifying how the stellar contraction leads to 
rotational spin-up and changes in interior structure for 
solar and intermediate mass stars, factors which com¬ 
bine to determine the initial conditions for the evolution 
of magnetic activity for the duration of the stellar life¬ 
time. Members of young stellar associations are usually 
dated by comparing their temperatures and luminosities 
to pre-main sequence tracks. However, errors and biases 
in placing stars on Hertzsprung-Russell (HR) diagrams 
lead to uncertainty in ages of different clusters, while 
large luminosity spreads among members of individual 
clusters enhance uncertainty r egard i ng age dating of any 
individual star (see review bv lSoderblom et aLll2014H . 

Further uncertainties in pre-main sequence age es¬ 
timates from HR diagrams are traced directly to the 
varied stellar loci predicted from different stellar evo¬ 
lution models and to the lack of consistency between 
predicted and observed stellar loci. These differences 
have been noted between low and intermediate mass 
stars (e.g. iHillenbrandl Il997t IHillenbrand et al.l 120081 
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(e.g. Ida Rio et al 1l20 lHlMurohv et alJl2013HKraus et al. ' 
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frequently have different estimates for age des pite an 


expectation that the y are c o-eval (lHartigan et al . 


Hartigan fe Kenvoi 


2003i : iKraus fe Hillenbrandl 


1994 : 


200ft 


Cennaro et all 12012 1. As a direct consequence, age es¬ 


timates depend on the considered range i n spectral type 
or stellar mass. IHillenbrand et ahl ( 2008 1 found that in 
many regions G stars appears ~ 2 — 5 times older than 
K stars, de pen din g on w hich pre-main sequence tracks 
were used. iPecaut et al.l (|2012fl quantified these differ¬ 
ences by measuring an age for the Upper Sco OB Asso¬ 
ciation of 11 ± 2 Myr, larger than the 5 Myr age mea¬ 
sured previously from pre-main sequence lo w mass stars 
(|Preibisch et all 120021 : iSlesnick et all 120061) . Similarly, 
in both younger and older clusters, analyses of main 
sequence turnoff and intermediate mass stars in color- 
magnitude diagrams yield ages that are twice as large as 
those obtained from low mass pre-main sequence st ars 
jM avne fc Navlorll2008l iNavloil1200ft iBell et aJll2013f ). 

Nearby stellar associations with ages of 5-40 Myr 
provide us with a laboratory to establish empirical 
isochrones, which can be used as a model-independent 
clock and as a benchmark for testing pre-main sequence 
models. Within a single association older than 10 Myr, 
the age spread expected from star formation processes is 
small compared to the absolute age. Member stars have 
traveled far from any remnant molecular material asso¬ 
ciated with their birth cloud and are unaffected by inter¬ 
stellar extinction. Most disks have already dispersed by 
these ages, so observational uncertainties related to ac¬ 
cretion and circumstellar extinction are negligible for the 
vast majority of stars in this age range. In recent years, 
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the growing number of proper motions and parallax dis¬ 
tances have led to a substantial improvement in the ac¬ 
counting of membership within you ng associations (e.g. 
iMalo et, aT1l2013t IKraus et, al.ll2014fl . The HR diagrams 
of these nearby moving groups should consequently have 
only a small age spread. 

At younger ages, placing stars on HR diagrams is 
more problematic because accretion and extinction af¬ 
fects temperature and luminosity measurements, and be¬ 
cause the affect of any genuine age spread on luminosities 
decreases (in dex) with age. Nevertheless, wide spreads 
in luminosities of members of young clusters have been 
interpreted as true age spreads. For example, in the ONC 
a 1.5 dex luminosity spread (Da Rio et al. 2010) led 
Reggiani et al. (2011) to claim an age spread of 1.5- 
3.5 Myr, although some of this spread may alternately 
be interprete d as a spread in radius from different accre - 
tion histories (iHartmann et al.ll997tlBaraffe et ahll2009f) . 
Improved spectral types and extinction measurements of 
two outliers in the ONC sample indicate that the lu¬ 
minosity spread may be overestimated dManara et al.l 
120131 ). If we consider ~ 2 Myr as an approximate age 
spread, then that would lead to 0.16 dex and 0.08 dex 
spreads in luminosities of 5 Myr old and 20 Myr old 
clusters of single stars (as assessed for 3500 K stars 
from the Baraffe tracks, see §2.5). These luminosity 
spreads are consistent with the 0.05 to 0.15 dex errors 
in placing young sta r s in the HR diagram ( see Table 1; 
iPreibisch et al.ll2012t iSoderblom et alil2014j l. 

In this paper, we use literature spectral types and 
photometry to construct empirical isochrones of lumi¬ 
nosity versus temperature of diskless low mass stars in 
nearby associations and compare them to pre-main se¬ 
quence isochrones. This approach is simil ar to the em¬ 
pirical iso c hrones derived by iHillenbrandl (| 20091 ) and by 
I Bell et all (| 2012 L 120131120141) , but is focused exclusively 
on low mass stars in nearby young associations and mov¬ 
ing groups and is motivated by a desire for deriving sim¬ 
ple formulations that approximate empirical isochrones 
in HR diagrams. In §2, we describe the literature data 
that is used to place stars on HR diagrams and the pre- 
main sequence evolutionary tracks used for age estimates. 
In §3 we develop empirical isochrones for a set of bench¬ 
mark young associations. In §4, the empirical isochrones 
are compared to isochrones predicted from models of pre- 
main sequence evolution to test evolutionary tracks and 
to obtain age estimates for the young clusters for each 
model. In §5, we discuss how the ages from different 
tracks compare with independent age measurements. 

2. SAMPLE, DATA, AND PRE-MAIN SEQUENCE 
EVOLUTION MODELS 

2 . 1 . Membership and Properties of Nearby Moving 
Groups 

The (3 Pic and Tuc-Hor Moving Groups and the e Cha, 
77 Cha, and TW Hya Associations are used in this paper 
as benchmark clusters to measure empirical isochrones of 
young regions. The results obtained from these bench¬ 
mark clusters are then applied to the Upper Sco OB As¬ 
sociation. Searches for membership in existing all sky 
surveys, along with analyses of proper motions and radial 
velocities, have substantially improved the confidence in 
membership for these associations. Meanwhile, the pro¬ 
liferation of accurate parallax measurements for the near¬ 


est groups has significantly reduced the uncertainty in lu¬ 
minosities arising from errors in distance measurements. 
When combined, these advances allow for a robust anal¬ 
ysis of the stellar loci of these clusters. 

In some analysis below, the two older associations 
(the /3 Pic and Tuc-Hor Moving Groups) and the three 
younger associations (TW Hya, 77 Cha, and e Cha) are 
grouped together. Gravity differences between these 
clusters are described in § 2 . 2 - 2 .3, and the similarity in 
ages amongst the older and younger groups will be es¬ 
tablished in §3. The binary census is discussed in this 
section, while the implementation of binaries is described 
in §2.4. 

Beta Pic Moving Group: Membership fo r the d Pic 
moving group is obtained from the analysis of lMalo et al.1 
(|2013L l2014bD . based on a > 90% probability calculated 
from proper motion, radial velocity, photometry, and 
parallaxes when available. Spectral types are obtained 
from their compilation of literature spectral types. Mem¬ 
bers with distances obtained from kinematic modelling 
are not used in this paper because many targets are close 
enough that the percentage uncertainties are too large 3 . 
Binarity was acco unted for based on the c ompilation of 
Malo ct al. (2013), which relied heavily on Janson ct al. 
( 20121 ) for VI-dwar fs. A few obj e cts an d binary informa¬ 
tion is added from lRiedel et al 'l (|2014T ) . 

Tuc-Hor: The membership and spectral t ypes for 
t he T uc-Hor Association is obtained from IKraus et all 
(120141) . The distances used here are from their kinematic 
analysis 4 . The spectral types used here are obtained from 
their spectral analysis rather than the SEDs. Binary in¬ 
formation is not included in our analysis of Tuc-Hor. 

e Cha Association: Membership, spectral types, and 
kine matic distances of e Cha members were obtained 
from iMurohv et al.l (120131) . This membership list is 
mostly consistent with b ut more complete than that of 
ILopez M arti et a ll (12013D. The ta rgets listed as mem- 
bers in Lopez Marti et all (120131) but excluded from 
IMurohv et al.l ( 20131) are~VW Cha and CM Cha, both of 
which have disks and would therefore not be included in 
our analysis, and 2MASS J12074597-7816064, which was 
identified by iFang et all (|2013l) as a proper motion out- 
lier. The kinematic distances from IMurphv et ahl (120131) 
are adopted here and have a median of 111 pc. The 
known bin aries are ac counted for following the descrip¬ 
tion in IMurohv et al.1 (120131) . though the binary census 
is poor. 

77 Cha Cluster: Membership in the 77 Cha clus- 
ter is based pr i mar ily on the initial identification by 
IMamaiek et al.l (119991) and on a follow-up census by 


3 The confirmed members (> 90% probability) with parallax 
distances occupy a tight stellar locus, with a luminosity scatter 
of 0.14 dex around the best fit line. Confirmed members between 
3300—5000 K with only kinematic distances have a luminosity scat¬ 
ter (standard deviation about the mean linear fit) of 0.36 dex (0.25 
dex if restricted to 3350—5000 K). 

4 Unlike the case for the /3 Pic Moving Group, the kinematic 
distances for Tuc-Hor produce a tight stellar locus, perhaps be¬ 
cause the average distance is larger so the fractional uncertainty 
is smaller. IKraus et al.l 12014 ) calculate that a 5% uncertainty in 
kinematic distance measurements applies to stars in their sample 
with a typical 5 mas yr~ 1 error in proper motion. T his claim is sup¬ 
ported by parallax measurements of five stars ((van Lccuwcnl l2007l : 
IRiedel e t al. 2014) that differ by 4% from the kinematic distances 
in Kraus et al. 
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Fig. 1.— The J-H and J-K colors for members of old (blue squares) and young (red circles, double circles indicate accretors) nearby 
associations, compared with colors from the BT-Settl models (log# = 4.0 is the solid line, log g = 5.0 is the dotted black line), the color 
from Pecaut & Mamajek (203, green line for young stars, dotted green line for main sequence stars), and the color adopted here (purple 
solid line for the older moving groups and purple dashed line for young associations). The J-H colors of (3 Pic Moving Group and Tuc-Hor 
Moving Group members are offset by ~ 0.05 mag from the colors of TW Hya, 77 Cha, and e Cha Association members. This measured 
gravity dependence is significantly less than predicted from the BT-Settl models for M-dwarfs. The J-H color is adopted by approximating 
the J-H color intermediate between the young and older clusters. 


iLuhman fe Steeghsl (|2004f) . The spectral types for t hose 
initial 18 members are obtained from ILuhman fc Steeghsl 
(12004) . Five dispers ed very low ma ss sta rs were found 
and characterized bv iMurphv et al.l (I2010D and are also 
included here. The a dditional members found by 
ILopez M arti et a l~ (120131) based on proper motions and 
the lCovino et, al.l ~ 19971) X-rav survey are excluded here 
because they are far from the cluster core and may there¬ 
fore be located at a different distance. The distance to all 
members is assumed to be 94 pc, the average Hipparcos 
parallax distance of p Cha and RS Cha (van Leeuwen 
2007; see discussion in Murphy et al. 2010). An ac- 
counting of spat i ally r esolved binaries was obtained by 
IBrandeker et al.l (|2006D and is applied here. 

TW Hya Association: Membership in the TWA 
is based on kinematic anal y ses dMalo et al.l 
2014bl: iWeinberger et al.l 120131 : iDucourant et al.l 


2013 


2014 


Riedel et al.l 1 201 4b Spectral types of most members 


a nd bina ri ty are obtained_ from the compilation in 
iHerczeg fc Hillenbrandl (12014 . The new equal luminos¬ 
ity binary TWA member SCR 1012-3124AB is also in¬ 
cluded here dRiedel et al.ll20l4 . Although parallax dis¬ 
tances are preferred, TWA stars with kinematic distance 
measurements are included in our analyses. For TWA 8A 
and 8B, a parallax distance distance of 44 pc is adopted 
here as the average of para llax distance s meas ured by 
IDucourant et al.l (12014 and iRiedel et al.l (12014 . Three 
members of the HD 98800 quadru ple system are included 
here with data from lLaskar et al.l (I2009D . despite the 
presence of weak IR emission from HD 98800 N. 

Upper Sco OB Association: The stellar census 
and spectral types for the Upper Sco OB Associ¬ 
ation are obtained from the compilation of mem¬ 
bers by ILuhman fc Mamai ekl (I2012T) . Multiplicity 
is correc ted based on several incomplete companion 
searches (IKraus e t al.l 120081: IKraus fc Hillenbrandl 120121 : 

iBiller et al.l 1201 ll lLafreniere et al.ll20l4 . The luminosi- 
ties estimated for the Upper Sco OB Association mem¬ 
bers are corrected for extinction (see §2.4). The aver¬ 
age extinction 5 of Ay = 0.46 magfor low mass stars 
and 0.58 mag for brown dwarfs is intermediate between 
the average extinction of Ay = 0.78 mag for A stars, 


Ay = 0.59 mag for G st ars, and Ay = 0.35 mag for F 
stars iPecaut et all ( 2012 1. Stars with Ay > 2 mag are 
atypical for Upper Sco and are not included in the HR 
diagram analysis. All stars in Upp er Sco are assu med to 
be located at a distance of 145 pc (Ide Zeeuwlll999l l. The 
cluster is 30 pc wide, and a depth of similar length likely 
contributes to the large l uminosity spread tha t has been 
measured for Upper Sco dSlesnick et al.ll2008 ). 

Stars without any WISE W2 or W3 photometry (2 
of 214 stars total between 3200-5000 K, 7 of 415 total 
brown dwarfs between 2800-3200 K) are included in the 
Upper Sco sample because their exclusion could bias any 
results at low masses against faint objects. 


2 .2. Colors of the Five Benchmark Clusters 

Photometry for members of the five benchmark 
clusters are obtained from the 2MASS JHKs survey 


(iSkrutskie et al 


2006), WISE Wl, W2, W3, W4 sur¬ 


vey dCutri et. all 2012f). the UCAC4 APASS all-sky BV- 
gri survey (iZacharias et al.l I2013IL a nd the GALEX 


NUV/FUV survey ( Bianchi et al.l I20TTT) . This collection 
of photometry is used in the following subsections to 
identify and exclude disks, to ensure that extinction is 
negligible (or correct luminosities for extinction, in the 
case of Upper Sco members), and to calculate bolomet- 
ric corrections (see also Appendix A). All objects are 
detected in JHK. Some objects lack photometry at other 
bands. 


TABLE 1 

Approximate Uncertainties in Stellar Properties 


Parameter 

T eff Range 

Uncertainty 

Temperature 

> 4100 K 

200 K or 2 SpT subclasses 

Temperature 

3750-4100 K 

125 K or 1 SpT subclass 

Temperature 

2800-3750 K 

75 K or 0.5 SpT subclasses 

Luminosity 

all 

~ 5% parallax distance 

Luminosity 

all 

10 - 20 % kinematic distance 

Luminosity 

all 

2-5% bolometric correction" 

“Or 0.02-0.05 

in mag 



5 The extinction for both low mass stars and brown dwarfs is 
0.66 mag if the photospheric J-K color is adopted from the older 
moving groups rather than the younger clusters. 
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Fig. 2.— The K-W3 color as a disk indicator. The samples 
analyzed here are restricted to diskless stars, identified by a K-W3 
color less than 0.5 mag above the best fit 2nd order polynomial 
(solid black line shows photospheric K-W3 color, dotted black line 
shows the cutoff for the presence of excess W3 emission). The red 
circles are stars in the three younger associations, the blue squares 
are stars in the two Moving Groups, and the small cyan circles are 
stars in Upper Sco. 

Broadband colors of objects in the younger TW Hya, 
r] Cha, and e Clia Associations are mostly consistent 
with those from the older /3 Pic and Tuc-Hor Moving 
Groups, which all form a tight stellar locus in infrared 
and optical colors (Figure[T]and Appendix). The contrac¬ 
tion to the main sequence corresponds with a gradual 
increase in gravity. Just as different gravities between 
mai n sequence and young stars yields color differ ences 
(e.g. lLuhman et aLll2003l : iPecaut fc Mamaiel3l2013h . cor¬ 
responding differences should be present between 5 and 
30 Myr old pre-main sequence stars. The J H color is 
redder for the younger associations relative to the older 
moving groups (Table [2] and right panel of Figure [T]), 
which is the result of gravity differences and possible mis¬ 
matches between SpT and temperature between the two 
samples. The J-K color is also ~ 0.03 mag redder for 
early-to-mid M stars in younger associations. Most em¬ 
pirical colors for the sample stars are also consistent with 
the in trinsic c olors for the 5-30 Myr stars as assessed by 
IPecaut fe Mamaieki (l2013lh 

The sample is restricted to diskless stars. Disks are 
identified by at least one of a K-W3 excess or K-W2 of 
0.5 mag above the photospheric level, as determined from 
the obseved stellar locus and from photospheric models 
(Figure O see also, e.g., Luhman & Mamajek 2012). A 
few stars may retain transition or debris disks that do not 
appear as a K-W2 or K-W3 excess, but these disks are 
unlikely to affect the J-band brightness of the star. The 
sources 2MASS J08014860-8058052, 2MASS J11550485- 
7919108, and 2MASS J02153328-5627175 fall below our 
disk criteria but are excluded in our analysis because of 
excess emission in the WISE W4 band or the detection 
of ac c retio n in optical spectra dMurphv et al.ll201ll 120131 : 
iCutri et al.ll2012f) . 

2.3. Bolometric Corrections and Temperature Scale 

The luminosities of stars in the five young, benchmark 
associations are calculated by applying bolo metric cor¬ 
rectio ns to 2MASS J-band magnitudes (Skrutski e et al.1 
120061 1. Any bulk offset between the actual and imple¬ 
mented bolometric correction would lead to an offset in 
luminosity and therefore age, while any temperature de¬ 


pendence in bolometric correction offset would change 
the relative huninosities/ages of stars in the same associ¬ 
ation but with different temperatures. These systematic 
uncertainties introduce problems when comparing differ¬ 
ent methods and when comparing observations to mod¬ 
els. In this subsection, we calculate bolometric correc¬ 
tions and subsequently assess temperature dependence 
in the uncertainties by comparing results from temper¬ 
ature and bolo metr ic c orrect ions derived here to those 
of IPecaut fc Mamajek! ( 2013 ). The Appendix describes 
the color corrections and constants that are used here to 
calculate fluxes and stellar luminosities. 

Our bolometric corrections are calculated from the BT- 
Settl model atmospheres (version CIFITS2011bc, which 
use a bundan ces fro m Caffau et al. 2011) with logg = 4.0 
llAllard et al.l[2012 f). Minor offsets relative to the J-band 
flux are applied to several spectral bands to account for 
differences between synthetic and observed photometry 
(see Figurefl3l). Spectral types ar e converted to tempera¬ 
t ure f rom the conversion scale of lHerczeg fc Hillenbrandl 
dMl, which were calculated by comparing low resolu¬ 
tion optical blue-red spectra to 2011 BT-Settl models. 
The model spectra and observed photometry are com¬ 
pared in magnitude space. The synthetic photometry 
is calculated by convolving the BT-Settl model spectra 
in photon units with filter transmission curves (see Ap¬ 
pendix A). The empirical colors are obtained from the 
photometry described in §2.2. The primary uncertain¬ 
ties in these empirical corrections are introduced by the 
lack of photometric coverage in the z and Y bands and 
in near-IR regions with poor telluric transmission, which 
are located between the zYJHK filters. 

Table [2] compares our spectral type-temperature con- 


TABLE 2 

Tempeartur.es and Bolometric Corrections versus Spectra 

Type 


SpT 

hhiT 

Tett 

PM13 C 

Here 

BC} 

PM13 C 

01d d 

J-H 

Young d 

PM13 C 

F5 

6600 

6420 

0.79 

0.85 

0.17 

0.17 

0.19 

F8 

6130 

6100 

0.95 

0.96 

0.22 

0.22 

0.23 

GO 

5930 

6050 

1.02 

0.98 

0.26 

0.26 

0.24 

G2 

5690 

5870 

1.10 

1.03 

0.30 

0.30 

0.27 

G5 

5430 

5500 

1.18 

1.16 

0.34 

0.34 

0.33 

G8 

5180 

5210 

1.26 

1.25 

0.40 

0.40 

0.39 

K0 

4870 

5030 

1.36 

1.30 

0.47 

0.47 

0.43 

K2 

4710 

4760 

1.41 

1.40 

0.50 

0.50 

0.49 

K5 

4210 

4140 

1.56 

1.58 

0.63 

0.63 

0.64 

K7 

4020 

3970 

1.62 

1.63 

0.66 

0.66 

0.66 

M0 

3900 

3770 

1.66 

1.69 

0.67 

0.68 

0.68 

Ml 

3720 

3630 

1.73 

1.74 

0.67 

0.70 

0.68 

M2 

3560 

3490 

1.78 

1.80 

0.65 

0.70 

0.67 

M3 

3410 

3360 

1.84 

1.84 

0.63 

0.70 

0.66 

M4 

3190 

3160 

1.93 

1.91 

0.58 

0.65 

0.62 

M5 

2980 

2880 

1.99 

2.01 

0.53 

0.59 

0.55 

M6 

2860 


2.03 


0.50 

0.55 


M7 

2770 


2.06 


0.48 

0.53 



a Following ^ecautj&JMamajek (2013), the zero point for the 
J-band magnitude is 3.129 x 10 —10 erg cm -2 s —1 A -1 and 
for the absolute magnitude is 3.055 x 10 35 erg cm -2 s —1 
b HH14 refers to Herczeg fc Hil lenbrand (2014) 

C PM13 refers to the 5-30 Myr stars in Table 6 
of ^ecaut^&^^amajek (2013) 

d J-H colors given for the median of diskless stars from the 
two old moving groups and three young associations 
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Fig. 3.— Predictions from different pre-main sequence evolutionary tracks. The top panels show how luminosity is converted to age for 
4200 K stars (roughly K5 spectral type, left panel) or 3400 K stars (roughly M3 spectral type, right panel). The bottom panels illustrate 
the temperature predictions at fixed luminosity. The effective temperatures are much higher in the DM97 models and, for L = L §, lower 
in the BCAH98 models. 


ve rsion a nd b olo metr ic correc tions to the 5-30 Myr stars 
of Pecaut_& - ^amajek ([2013). In general, temperatures 
in Pecaut fc Mamajekj (12013H are 100-200 K cooler than 


those presented here. These errors are roughly assessed 
from the uncertainty in assignin g spectra l ty pes fr om low 
resolution optical spectra. iPecaut fc Mamaiekl (|2013T ) 
also include K8 and K9 spectral types. Their SpT- 
temperature conversion extend only to M5, so results 
for br own d warfs ( §3.2 ) are limited to the conversions of 
iHerczeg fo Hillenbrandl (12014H . When evaluated at the 
same temperature, the bolometric corrections measured 
di rect ly from the BT-Settl models with thejCaffan^tjjl 


(2011) abundances reproduce the IPecaut fe Mamaiek 


(20131) bolometric corre ctions to within 1 % at all spec¬ 


tral types, even though IPecaut fc Mam aiek 


models with abundances from lAsphmd et at 


ferences between our bolometric corrections and those 
of Pecaut et al. are introduced by differences in the 
SpT-temperature conversion and by our tweaking of the 
BT-Settl models to reproduce the measured colors, rel¬ 
ative to J. The bolometric correction versus spectral 
type is also different because of the different spectral 
type-temperature conversion scale. As a consequence, 
when compared by temperature, bolometric corrections 
are similar for M-dwarfs and ~ 2% smaller here for K- 
dwarfs. The practical result of these differences is that 
temperatures are ~ 50 — 150 K higher and luminosi¬ 
ties are ~ 1.5 — 4.5% b righter relative to those calcu¬ 
lated from lPecaut fe Mamaiekl (1201311 bolometric correc¬ 
tion and temperature scales. 


The analysis of stellar temperatures and luminosi- 
t ies f or 320 0 -500 0 K sta rs (§3.1) is repeated for the 
IPecaut fc Mamaiek (12013[ 1 SpT-tem perat ure conversions 
and bolomet ric correcti ons. Both IPecaut fc Mamaiekl 
(|20H and IHerczeg &; Hillenbrandl (120141) base their 
SpT-tem perature conversions on the BT-Settl models, 
although lTottle fc Mohantvl (|2015H find discrepancies be¬ 
tween these models and near-IR colors at both early and 
late-M spectral types. While the SpT-temperature con¬ 
version is still a significant source of uncertainty, the 
main results of this paper are robust to systematic un¬ 
certainties in these scales. 

2.4. HR Diagram Methodology 

Stellar luminosities are calculated from 2MASS J-band 
brightness, literature spectral types, and stellar distances 
calculated from parallax measurements or from asso¬ 
ciation kinematics. The J-band is used for luminos¬ 
ity calculations because of the uniform reliability of the 
2MASS survey, and because the J-band is near the peak 
of the spectral energy distribution for late type stars, 
which minimizes the uncertainty in bolometric correc¬ 
tion caused by spectral type errors. The stars are disk¬ 
less, so accretion and disk extinction do not affect the 
observed J-band photometry. Stars in each association 
are assumed to be co-eval. 

Table [T] lists the uncertainties applied here for luminos¬ 
ity and temperature estimates. Luminosity uncertainties 
are calculated from an assumed 5% uncertainty from a 
combination of photometric and bolometric correction 
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Fig. 4.— A comparison of the new Feiden (red) and BHAC- 
2015 (blue) tracks with an improved treatment of convection. The 
solid lines show isochrones at log age/yr=6.0, 6.5, 7.0, 7.5, and 8.0, 
with 6.0, 7.0, and 8.0 in bold. The vertical lines show the early 
evolution of 0.12, 0.3, and 1.0 M© stars from log age/yr=5.7 to 
8.0. Although the physics implemented into these modern tracks 
are different, they produce evolutionary tracks and isochrones that 
are extraordinarily similar. 

errors and from distance uncertainties, typically < 5% for 
parallax distances (5% distance uncertainty is assessed 
for all 77 Cha Association members), 5% for kinematic 
distances for Tuc-Hor members (see footnote 3), and 10% 
for kinematic distances of other regions. Bolometric cor¬ 
rections for the J-band are described in detail in the next 
section and are more uncertain for cooler stars with large 
opacities in molecules. The depth of Upper Sco adds 
uncertainty in the distance and therefore luminosity for 
individual members, which contributes to its luminos¬ 
ity spread. Temperature uncertainties are based on esti¬ 
mated errors and methodological differences in spectral 
typing (see discussion in Herczeg & Hillenbrand 2014). 
Systematic differences in spectral typing between differ¬ 
ent studies will also affect the HR diagrams constructed 
in this paper. 

In multiple systems, all members with resolved spectral 
types are included in this analysis. Unresolved spectral 
types are included if the stars are spatially resolved in J, 
H, or K-bands. In these cases, the primary is assumed 
to have the system spectral type. The J-band magni¬ 
tude of the primary is calculated by assuming that the 
brightness difference is constant in the near-IR (the J-H 
and J-K colors roughly constant for M-dwarfs). Opti¬ 
cally resolved or double-lined spectroscopic binaries are 
included if both stars are similar brightness. Single-lined 
binaries are excluded. 

Extinction is negligible and assumed to be negligi¬ 
ble for all diskless stars in the five benchmark associ¬ 
ations (see colors in §2.2 and the Appendix). Stars in 
the younger associations have a slightly bluer g-J color, 
which is the opposite of expectations if extinction were 
important and supports the assumption that extinction 
is negligible. The four stars i n e C ha wi th Ay >1.5 mag, 
as calculated by iMurphv et al.i (120131) . all have disks, 
which may either produce extinction or may alter the 
near-IR colors sufficiently to have led to extinction in 
their calculations. 

However, extinction affects luminosities of stars in the 
Upper Sco OB Association. A color excess E(J-K) is 
calculated for each star here from by measuring the dif¬ 
ference between the observed J-K color the J-K color 


adopted for the young associations, as shown as the 
dashed purple line in Figure [l] The extinction Aj is 
then calculated and corrected for by assuming a total- 
to-selective extinc tion of Aj = 0.2827W , and Aj — 
1.66 E(J—K) from lRieke fc Lebofskvl (|1985f) . Ifthetotal- 
to-selective extinction Ry = 4 instead of 3.1, the Aj 
values would increase by 0.0 7%, as calcula ted from the 
extinction curves of lWeingartner fe Drainel (120011 ). 

Young K and M stars without disks have extinctions 
accurate to ~ 0.2 mag, although the precise error in ex¬ 
tinction depends on the error in spectral type. Stars with 
extinctions Ay > 2 mag are excluded from our analysis. 
Negative values for extinction are retained and applied 
for statistical purposes. 


2.5. Models of pre-main sequence evolution 

The empirical isochrones derived in §3 are compared 
topre-main sequence isochrones in §4. iStassun et al.l 
( 2014 ) summarizes currently available pre-main sequence 
evolutionary tracks. The evolutionary models used here, 
D’Antona & Mazzitelli (1997, hereafter DM97), Siess et 
al. (2000, hereafter Siess), Baraffe et al. (1998, 2003, 
hereafter BCAH98), Dotter et al. (2008, hereafter Dart¬ 
mouth), Tognelli et al. (2011, hereafter Pisa), and PAR- 
SECvl.2s (Bressan et al. 2012, Chen et al. 2014; her- 
after PARSEC, also frequently referred to as Padova), 
are selected based on their availability and prevalence in 
the pre-main sequence literature. The BCAH98 models 
have mixing lengths a = 1.0 (BCAH98) or mixing length 
a = 1.9 (hereafter BCAH98-1.9). The PARSEC models 
are selected for their artificial tweaks to the T — r re¬ 
lationship that reproduce the color-magnitude diagrams 
of very low mass stars in Praesepe and M67. For stars 
below 4000 K, stars of the same age and mass are much 
cooler and are only slightly less luminous (because of 
larger radii) in the PA RSECv l .2s models r elative to the 
PARSECvl.l models ( Bressan et al.ll2012f ). The PAR- 
SECvl.l models are used only in Tableland are not 
discussed further. 

In addition to these sets of models, we include anal¬ 
ysis of newly available tracks calculated by Baraffe et 
al. (2015, hereafter BHAC-2015) using the BT-Settl at¬ 
mospheric models and a prescription of convection based 
on 2D/3D models, which is an improvement from the 
mixing length prescription for convection. We also in¬ 
clude new Dartmouth tracks, which improve the thermal 
structure in the atmosphere by calculating fitting to large 
optical depths (Feiden et al. in prep., private communi¬ 
cation, hereafter F e iden) . The Feiden models use the 
lOrevesse fc Sauval (1 1 9981 ) a bundances with PHOEN IX 
AMES-COND atmos pheres (lHauschildt et al .11 1 999al fbl) . 
IStassun et al.l (|2014l) includ e prelimina r y descr iptions 
of the Feiden tracks, while iMalo et~akl (l2014al) use a 
magnetic version of the Feiden tracks. The BHAC- 
2015 and Feiden models predict very similar tempera¬ 
ture/luminosity evolution for low mass stars and brown 
dwarfs, which then yields very similar isochrones (Fig¬ 
ure 13}. 

All model tracks use solar abundances 6 , which is con¬ 
sistent with abundance mea surements of y oung stars in 
the solar neighborhood (e.g. iD’Orazi et aDl201lD . Only 
the BCAH98 and the DM97 pre-main sequence tracks 
are available for masses <0.1 M 0 , while the Pisa mod¬ 
els are only available for masses > 0.2 M 0 . The BCAH98 
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Fig. 5.— The evolution of temperature and luminosity for a 0.4 Mq star, as predicted by different pre-main sequence models. The 
temperature predictions are disparate at the 10-15% level among models at all ages. However, the luminosity evolution are similar after 2 
Myr for all models of the same mass, despite the different effective temperatures. The BHAC-2015 tracks follow the Feiden tracks and are 
not shown here. 


models with a = 1.0 are available for stars less massive 
than < 1.4 Mq while only models between 0.6 — 1.4 Mq 
have been calculated with a = 1.9. The PARSEC models 
are only calculated to 0.1 Mq, but extend to low tem¬ 
peratures (2800 K at 1 Myr, 2500 K at 10 Myr) because 
the tweaks to the T — r relationship reduce the effective 
temperature of very low mass stars. 

Figure [3] shows the conversion to age at fixed lumi¬ 
nosity and at fixed temperature for the different pre- 
main sequence evolutionary models. The DM97 m odel 
track s are significant outliers (see also, e.g . [Hillcnbrand 
Il997t iBell et al.ll2012t lAndrews et aI1l2013 ). Their pho- 
tospheric temperatures are hotter for stars of the same 
mass, so the corresponding mass of K and M stars 
is much lower than obtained from other evolutionary 
tracks. The measured luminosity for a star of that mass is 
therefore high, leading to a young age. At ages < 3 Myr, 
the Dartmouth tracks also yield ages that are younger 
than the other tracks. In contrast, the BCAH98 tracks 
assessed for 4200 K stars are outliers in yielding older 
ages. 

All models predict similar luminosities for stars of the 
same mass for ages older than 2 Myr, despite differ¬ 
ent temperatures (Fig. [5]). With accurate age dating of 
young clusters, masses may be inferred from the lumi¬ 
nosity alone and with consistent results between different 
model tracks. 

3. EMPIRICAL ISOCHRONES FOR THE NEARBY YOUNG 
MOVING GROUPS 

With an accurate theory of pre-main sequence evolu¬ 
tion, the estimated luminosity of a star at a given es¬ 
timated temperature could be directly converted to an 
age. However, scatter in L and T e g among young cluster 
members provokes concerns regarding the reliability of 
age estimates for single stars. Moreover, differences in 
the slope between observed and model isochrones lead to 
ages that depend on spectral type. Our approach here 
focuses on comparing the empirical stellar locus of dif¬ 
ferent associations and interpreting the offsets as relative 
age differences. Using the stellar loci constructed for 
young associations, we first focus on empirical trends in 

6 Changes in solar abundance measurements means that these 
models may use different abundances and are therefore not directly 
comparable. 


log L versus T e g (this section) before conducting detailed 
comparisons to models (§4). 

Figure [6] shows the combined HR diagrams for the five 
benchmark young associations. The stellar locus between 
3500-5000 K has a similar shape for all clusters, with 
offsets in luminosity and temperature. Stars in the older 
associations have systematically fainter luminosities and 
warmer temperatures relative to young associations. The 
slope of the locus is steeper for very low mass stars and 
brown dwarfs 7 (T < 3200 K for the young clusters and 
< 3400 K for older clusters) relative to the locus for solar 
mass stars. The location where the break in slope of 
the temperature-luminosity isochrone occurs is referred 
to here as the inflection point in the isochrone. 

The observed HR diagrams are consistent with the ba¬ 
sic shape predicted from isochrones obtained from pre- 
main sequence evolutionary models (see comparisons in 
§4 and Fig. ED- The gradual decrease in luminosity from 
5000 to ~ 3300 K is roughly consistent with all pre-main 
sequence evolutionary models. At temperatures cooler 
than ~ 3300 K, the steep dropoff in luminosity is also 
roughly consistent with predictions from tracks that cal¬ 
culate pre-main sequence evolution for stars < 0.2 Mq. 
Based on the HR diagrams, the inflection point where 
the slope steepens occurs at ~ 3400 K for the older 
two associations and ~ 3200 K for younger associations. 
The temperature intervals are therefore different in fits 
to stellar and brown dwarf loci of the older and younger 
associations. 

The similarity in stellar loci suggests that the mem¬ 
bers may be analyzed together after accounting for age- 
dependent differences between the clusters. In the follow¬ 
ing subsections, we combine the loci from different associ¬ 
ations into a single stellar locus to calculate an empirical 
isochrone, which can then be generically applied to any 
young cluster. Combining data from different clusters 
yields results that are more robust than those calculated 
for individual, sparsely populated associations. When 
comparing cluster ages, scatter in data points and differ¬ 
ent spectral type distributions of members can introduce 
systematic errors in relative ages. 

7 These objects include very low mass stars and brown dwarfs 
and are termed brown dwarfs in this paper for ease of terminology. 
Object masses are not calculated in this paper, and the distinction 
between stars and brown dwarfs is not relevant at young ages. 
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Fig. 6. — HR diagrams for nearby young associations. The top diagram shows the luminosity and temperature for members of the 
benchmark five associations (green diamonds for the Tuc Hor Moving Group, red circles for the j3 Pic Moving Group, black asterisks for 
the TW Hya Association, blue squares for the i] Cha Association, purple triangles for the e Cha Association, and small cyan circles for the 
Upper Sco OB Association). The bottom panels show the same figure, with the stars offset in luminosity (left) to calculate a best fit to 
the 3200-5000 K stars (black solid line) and offset in temperature (right) to calculate a best fit to the 2800-3200 K brown dwarfs (black 
solid line, extends to < 2800 K because some temperatures are shifted into this range). 


Once a generalized fit to temperature and luminosity 
is obtained, scaling the slopes to each cluster in the lu¬ 
minosity direction then provide self-consistent compar¬ 
isons of relative ages that furthermore are robust to cer¬ 
tain types of scatter in the temperature-luminosity data 
points. For example, linear fits to the TWA and Tuc- 
Hor populations yield median luminosities of log L/L@ = 
—0.34 and —0.70, respectively, for 4200 K stars, or —1.00 
and —1.34 for 3400 K stars 8 . These luminosities can be 
compared among clusters to infer relative ages in log L 
as a fuction of log T. 

8 These luminosities differ from the final values listed in Table [4] 
because these values are calculated from best-fit slopes for the two 


The loci of stars above and below the inflection point 
(3200 K for young regions, 3400 K for the older regions) 
are discussed separately in §3.1 and §3.2. An underlying 
assumption to this approach is that the slopes of stellar 
loci are constant longward and shortward of the inflec¬ 
tion point and that the slopes do not change with cluster 
age. These sections include a discussion of whether this 
simplifying assumption is correct. For brown dwarfs in 
particular, the logL — T slopes are different for well- 
populated loci of Upper Sco and Tuc-Hor. The location 
of the inflection point ranges from 3000-3500 K and is 

clusters, while the final values are obtained from scaling a line with 
a slope that best fits all regions together. 
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TABLE 3 

Slopes of Linear Fits to Young Stellar Isochrones 




3200-5000 K 

2800-3200 K 

Empirical Fit 


Slope 

StD 

Slope 

StD 

Combined 


0.76 ±0.07 

0.15 

2.05 ±0.15 

0.18 

Young Groups Only a 

0.76 ± 0.06 

0.16 

3.1 ±0.8 

0.20 

Old Groups Only*’ 

0.85 ± 0.08 

0.14 

2.05 ±0.18 

0.18 

Tuc Hor 


0.83 ±0.12 

0.14 

2.28 ±0.25 

0.18 

BPMG 


0.91 ± 0.12 

0.14 

1.87 ±0.38 

0.16 

TWA 


0.82 ± 0.15 

0.17 

d 


e Cha 


0.73 ± 0.08 

0.17 

d 


77 Cha 


0.93 ± 0.27 

0.06 

d 


Upper Sco 


0.85 ± 0.07 

0.18 

3.5 ± 0.8 

0.27 

Model 

Age c 

3400-5000 K 


2800-3200 K 


BCAH98 

3 

1.04 


1.97 


BCAH98 

20 

0.89 


1.62 


BCAH98-1.9 

3 

0.69 


- 


BCAH98-1.9 

20 

0.76 


- 


BHAC-2015 

3 

0.68 


1.49 


BHAC-2015 

20 

0.72 


1.84 


Pisa 

3 

0.71 


- 


Pisa 

20 

0.76 


- 


Dartmouth 

3 

0.73 


- 


Dartmouth 

20 

0.74 


- 


Feiden 

3 

0.75 


1.57 


Feiden 

20 

0.84 


2.02 


Siess 

3 

0.87 


- 


Siess 

20 

1.00 


- 


DM97 

3 

0.70 


2.46 


DM97 

20 

0.85 


2.94 


PARSEC 

3 

0.62 


1.15 


PARSEC 

20 

0.51 


1.32 


PARSECvl.l 

3 

0.83 


- 


PARSECvl.l 

20 

1.09 



- 


Slope (in log L per 10 a K), uncertainty, and standard deviation 

of fits to individual, combined, and model isochrones 
a ri Cha, e Cha, and TW Hya Associations 


b /3 Pic and Tuc Hor Moving Groups 
c Age of the model isochrone in Myr 
d Insufficient number of data points 

discussed in §3.3. 

3.1. Empirical Isochrones for low mass stars 

In this subsection, we describe global fits to the 3200- 
5000 K region (3400-5000 K for the /3 Pic and Tuc-Hor 
Moving Groups), which corresponds to M3-K0 spectral 
types. Since the empirical loci in this temperature range 
are similar, a global linear fit to the entire dataset pro¬ 
vides a consistent assessment of the relative luminosity 
evolution of a given cluster. In fits to the combined 
HR diagram of the five benchmark associations, the free 
parameters are the slope of the line and a different y- 
intercept for each cluster. For illustrative purposes, the 
different y-intercepts are described here as luminosity off¬ 
sets. 

The lower left panel of Figure [ 6 ] shows the HR diagram, 
with each association shifted by the respective difference 
in y-intercept to match the location of the TWA. The em¬ 
pirical isochrone for the combined HR diagram of nearby 
associations has a best fit slope of 

log L/Lq oc (0.76 ± 0.07) T/10 3 K (T infl - 5000 K), 

. 

based on linear regressions that includes errors in T 
and L (as describ ed in §2.3) and calcu lated using mp- 
fitexy.pro in IDL dWilliams et al.l P20 1 01) . The inflection 


point Tj n fl used for the fits here is 3200 K for the younger 
groups and 3400 K for the older groups. The 0.07 dex un¬ 
certainty includes a 0.05 dex statistical uncertainty and 
a ~ 0.05 dex uncertainty based on different linear fit¬ 
ting techniques, how the data points are weighted, and 
the inclusion or exclusion of deviant data points. This 
uncertainty does not include possible errors in the spec¬ 
tral type-temperature conversion or in bolometric cor¬ 
rections, which may depend on spectral type. The stan¬ 
dard deviation in luminosity between the observations 
and best linear fits is ~ 0.15 dex. 

The same an alysis calculated from the 
[Pecaut fe Mamaieki (12013H spectral type-temperature 
conversion and bolometric corrections produces a slope 
of 0.72±0.04 log L per 10 3 K, which is formall y consistent 
with the results when using the iHerczeer fo Hillenbrandl 
1 2014 F1 conversions. The different slopes would yield a 
change in logT of 0.07 dex between the hottest and 
coolest stars in the applicable temperature range, which 
is smaller than other systematic errors. 

Table m lists the best fit slopes for linear fits to the 
the combined data set, to individual clusters, and to the 
clusters combined into the young and old groups. When 
shifted based on the best fit offsets, the three youngest 
benchmark associations have steeper slopes than the 
older clusters, which suggests some evolution with age. 
The difference in slopes between the young and old clus¬ 
ters, 0.76 and 0.85 logL per 10 3 K, may lead to system¬ 
atic differences as large as 0.14 dex for stars at 3400 K 
compared with 5000 K. 

3.2. Empirical Isochrones for Young Brown Dwarfs 

In this subsection, we describe global fits to objects 
cooler than the inflection point in the isochrones (2800- 
3400 K for the Tuc-Hor and /3 Pic Moving Groups, 2800- 
3200 K for the TW Hya, 77 Cha, and e Cha Associations). 
For illustrative purposes, the different y-intercepts are 
visualized as temperature offsets so that the loci over¬ 
lap (lower right panel of Fig. [ 6 ]). This visualization is 
practical because the inflection point occurs at different 
temperatures but approximately the same luminosity (to 
within a few tenths of dex; see also similar behavior in 
the model isochrones in Figs. IH and fTTl) . 

The total combined fit to our benchmark clusters yields 
a slope of 

log L/Lq oc (2.05T0.15) T/10 3 K (T infl -5000 K), (2) 

where the temperature range goes from 2800-3200 K for 
the younger clusters and 2800-3400 K for the older clus¬ 
ters (as also described above with Tj n fl). The error in 
the slope is calculated as in §3.1. The fits are dominated 
by members of the Tuc Hor and /3 Pic Moving Groups 
because those regions have more identified brown dwarf 
members. Fits are calculated for the combined locus of 
young associations but are not calculated for each indi¬ 
vidually because each association has only a few data 
points. Table[3]lists slopes for the individual regions and 
the combined HR diagram. 

The coolest known objects in these regions (< 2800 K) 
are fainter than expected from these linear fits. They 
are excluded in this analysis, but their inclusion would 
significantly steepen the slopes of this lines. This dis¬ 
crepancy suggests that these fits may not be robust to 
the selection of the low temperature cutoff for this fit. 
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Fig. 7.— The HR diagram for young (red circles) and old clusters (blue squares), with members of the j3 Pic Moving Group shifted to 
overlap with the Tuc Hor stellar locus, while 77 Cha and TW Hya Association members shifted to overlap with the e Cha Association stellar 
locus. The blue and red lines show the best fits to the combined low-mass star and brown dwarf locus, also shifted to match the young 
and old groups. 


The measured uncertainty in the linear fit of 0.15 
log L/Lq per 10 3 K leads to a 0.06 dex uncertainty in 
relative luminosity between a 3200 K and 2800 K brown 
dwarf. The difference in slopes between the older and 
younger associations likely leads to 0.09 dex uncertain¬ 
ties over the same regions. This difference is likely real 
evolution in the brown dwarf loci and not related to un¬ 
certainty in the fits. In §3.4, comparisons to the Upper 
Sco OB Association demonstrate that the slope becomes 
shallower with age. The results from linear fits in pro¬ 
ducing a luminosity for a 3000 K brown dwarf may still 
be robust as long as the known membership is well dis¬ 
tributed between 2800-3200 K. However, this relation¬ 
ship should only be used with caution. 

Th e fit to t h e brow n dwa rf locus is not recalculated 
for iPecaut fc Mamaiek l (j2013h parameters because their 
bolometric corrections are not calculated for spectral 
types later than M5. 

3.3. Combined Isochrones between 2800-5000 K 

Subsections §3.1 and 3.2 describe empirical fits to HR 
diagrams, as illustrated in Figure [6] The regions warmer 
and cooler than the inflection point in the temperature- 
luminosity locus were fit separately, with one slope but 
different y-intercepts for each association. The shifts ob¬ 
tained from the y-intercepts are visualized as luminosity 
shifts, while the fits to the cooler stars are visualized as 
temperature shifts. 

The location of the inflection point, which motivates 
the separate fits to the two temperature ranges, depends 


on the cluster age. In this subsection, we compare the 
independent results from the fits to low mass star and 
brown dwarf isochrones to describe how the empirical 
lines combine to describe the low mass locus of young 
clusters. For illustrative purposes only, the combined 
isochrones are plotted against the young and old associ¬ 
ations in Figure [7] 

Figure [5] shows the luminosity of 3000 K brown dwarfs 
versus 4200 K stars from our independent fits to the 
young clusters with a best fit line of 

\og(Li 2 oo/L q) = —1.02+ 1.54k>g(L 3 ooo/.£©)- (3) 

Combining the line of equation 3 with the two slopes re¬ 
ported in equations 1 and 2 fully describes the evolution 
of the empirical locus of single low mass stars and brown 
dwarfs in clusters. The isochrone for low mass stars falls 
in luminosity more quickly than that for brown dwarfs. 
The inflection point where the temperature-luminosity 
slope changes is determined by the intersection of the 
stellar and brown dwarf locus described by this line. 

As a further test of our method, Figure HI also shows 
preliminary results of its application to younger clusters 
with proportionally higher luminosities. The data are ob- 


tained for IC 348, Chameleon I, Lu 

pus I. Taurus, a Ori. 

A Ori, and the ONC ( 

jiihman et al. 

2003; Lulnnan 2007t 

Lulnnan et al. 2010; 

Rebull et all 

2010: Rieliaco et al. 

2011; Bayo et alj 2011]; Ida Rio et al.1 2012(). Although 


the errors in placement of individual stars on the HR 
diagram are larger for younger stars with more influence 
from disks and heavier extinction, their mean huninosi- 
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log L/L 0 (4200 K) 

Fig. 8.— The relationship between the luminosity for 4200 K and 
3000 K stars based on the lines fit to the individual clusters. The 
blue and red circles are the old and young associations (including 
Upper Sco), respectively, and are best fit with the solid black line. 
The purple squares are preliminary results from younger clusters, 
most of which are fit by an extrapolation of the best fit line to 
higher luminosities. 

ties are encouragement that our method could be applied 
to younger regions in order to constrain the luminosity 
evolution for comparison to theoretical models. 

3.4. Application to the Upper Sco OB Association 

Figure[9]shows how the fits derived above are applied to 
the HR diagram of the Upper Sco OB Association. Lines 
with slopes that best fit the loci of low mass stars (3200- 
5000 K) and brown dwarf/very low mass stars (2800- 
3200 K) are scaled to the measured isochrone. 

The initial sample includes stars with K-W3 color con¬ 
sistent with diskless stars and also stars without WISE 
detections. An initial scaling to the 3200-5000 K region 
(186 stars) has logl^oo/T© = —0.28 with a standard 
deviation of 0.26 dex, much larger than that obtained 
for the other clusters. Several stars are severe outliers, 
defined here as 0.5 dex from the logL calculated in the 
initial fit. Our adopted fit of 176 stars excludes these se¬ 
vere outliers. We calculate log L 4200 /T© = —0.23 with a 
luminosity scatter of 0.18 dex. Fitting stars from 3400- 
5000 K yields the same log £ 4200 - This luminosity spread 
of 0.18 dex is slightly larger than that calculated for the 
other associations, with extra scatter that may be caused 
by a real age spread in the cluster or, as suggested by 
[Slesnick et al. (200 8(1. a ~ 30 pc depth of the region. 

The brown dwarf locus of Upper Sco at 2800-3200 K 
has a much steeper slope than the linear fit to the 5-40 
Myr old sample, which points to problems in our sim¬ 
plifying assumption that the slope of the 2800-3200 K 
object locus does not evolve on these timescales. The fit 
to the slope of combined loci of the five benchmark clus¬ 
ters is dominated by Tuc-Hor members. However, the 
slopes from the younger associations are steeper than 
for the older associations. These differences suggest that 
the empirical brown dwarf locus is better described by 
a swinging gate rather than a single slope. This same 
swinging gate behavior is also suggested by the theo¬ 
retical isochrones covering this low temperature regime 
(Figure 9), where the predicted slopes appear to flatten 
with increasing age. In particular, the brown dwarf loci 
is remarkably well described by the DM97 isochrones and 
is also consistent with the BHAC-15 isochrones. 

While a more in-depth analysis is beyond the scope of 


TABLE 4 

Luminosity Ordering of Young Clusters 


Stellar 

Ass’n 

d c 

Low mass stars 0 

N log L4200 StD 

N 

Brown Dwarfs^ 
log L3000 StD 

Upper Sco 

145 

176 

-0.25 

0.18 

286 

-1.51 

0.27 

e Cha 

111 

14 

-0.28 

0.17 

4 

-1.42 

0.23 

77 Cha 

94 

7 

-0.32 

0.06 

5 

-1.61 

0.16 

TW Hya 

(57) 

22 

-0.40 

0.17 

3 

-1.51 

0.22 

P Pic 

(30) 

16 

-0.65 

0.14 

11 

-1.91 

0.15 

Tuc Hor 

(51) 

48 

-0.72 

0.14 

89 

-2.20 

0.18 


“From linear fits to 3200/3400-5000 K stars, see §3.1 
6 From temperture shifts for stars of 2800-3200 K, see §3.2 
“Distances in () are median association distances 


the simple-minded spirit of this paper, these problems 
with the brown dwarf isochrones lead us to focus sub¬ 
sequent analysis of the correspondence between stellar 
mean stellar luminosities and stellar ages on the 3200- 
5000 K stellar locus. 

3.5. Ordering Clusters by Luminosity 

In this section, we apply the analytical approximations 
of isochrones from §3.1-3.2 to the set of young clusters 
to order them in luminosity, as a proxy for age. The fits 
to luminosities for 3200-5000 K stars (§3.1) are likely 
more accurate than for the 2800-3200 K brown dwarfs 
(§3.2) because an analysis of the brown dwarf locus relies 
heavily on the two older moving groups and may not ac¬ 
curately describe the brown dwarf locus of younger clus¬ 
ters. 

Table [I] shows the ordering of clusters by luminosity 
for the 3200-5000 K fits. The five benchmark associa¬ 
tions are ordered in luminosity as expected from previ¬ 
ous studies (e.g., see discussions in Mamajek 2008; Hil- 
lenbrand 2009; Kraus et al. 2011 and Fang et al. 2013). 
The Upper Sco, e Cha and p Cha Associations are the 
youngest, followed by the TW Hya Association, the /3 Pic 
Moving Gro u p, and the Tuc-Hor Moving Group. Use of 
iPecaut fo Mamaiekl (2 013! ) temperatures and bolometric 
corrections would increase association logi^oo by 0.02- 
0.03 dex. 

The formal uncertainties in the luminosity (age) or¬ 
dering are calculated from a combination of (1) the stan¬ 
dard deviation divided by the square root of the number 
of stars in the sample, and (2) systematic uncertainty 
in distances. Methodological differences in luminosity 
calculations and in the census of multiplicity also can af¬ 
fect the luminosity ordering of clusters. For the specific 
benchmark clusters analyzed here, the age ordering is 
secure, with the possible exception of the rj Cha Associ¬ 
ation. To swap the TWA and g Cha Associations would 
require a systematic change in L4200 by 0.08 dex (20%). 
For young clusters, the primary culprits to explain such 
a luminosity change of 0.1 dex would be systematic er¬ 
rors in the distance (12% change), extinction (change of 
Ay = 0.9 mag), and sp ectral type (change of 140 K). 

I Weinberger et al.l (j2013l ) argue that TWA 9A and 9B 
should be excluded from the TWA census because they 
are much fainter than expected, which would increase the 
cluster luminosity by 0.04 dex. These comparisons are 
also affected by differences in binary accounting and by 
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Fig. 9.— The HR diagram of the Upper Sco OB Association, with empirical linear fits scaled in two temperature ranges to the stellar 
locus. Data points used for fits are solid circles. Objects with that are cooler than 2800 K or are more than 0.8 dex fainter than the median 
are not fit (empty circles). Objects with disks (cyan plus signs) are also excluded in the fits. The black lines show isochrones of log age 
= 6.0, 6.5, 7.0, 7.5, and 8.0, with 6.0, 7.0, and 8.0 in bold. The purple line shows the empirical isochrone for low mass stars and brown 
dwarfs, scaled to the Upper Sco OB Association data. 



Fig. 10.— The measured slopes of luminosity versus tempera¬ 
ture for each association (blue circles for 3200-5000 K stars with 
luminosities assessed at 4200 K; red circles for brown dwarfs with 
luminosities assessed at 3000 K) compared with slopes predicted 
from Feiden, BHAC-15, Pisa, DM97, PARSEC, and Siess tracks 
evaluated at 20 Myr (see Tables 3-4 for values). The solid lines 
show the luminosity interval over which the slopes apply for both 
the young and old associations in this paper. The dotted lines 
show the intermediate regime where the inflection point occurs. 
The brown dwarf points include one slope and luminosity calcu¬ 
lated for the combined young associations. 

accretion, if samples include disk sources. 

4. COMPARING EMPIRICAL AND MODEL ISOCHRONES 

The relationship between stellar luminosity and effec¬ 
tive temperature has been measured for six young asso¬ 
ciations. Because these stars are contracting pre-main 
sequence stars, the radius and therefore luminosity de¬ 
creases with time. The luminosity of stars at a fixed 
temperature is therefore used in §3.5 as a proxy for age. 
The conversion from luminosity to age depends on the 
adopted choice among the many evolutionary models, as 
shown in Fig. 3, and on the stellar temperature. At a 
fixed temperature, the age estimates differ between the 
many evolutionary models. Mismatches between empiri¬ 
cal and model isochronal slopes lead to different age esti¬ 
mates along an empirical isochrone when a given model 
family is applied to stars with different temperatures. 

Figures [lOlfTlI and Table [3] compare the empirical 
isochrones to the pre-main sequence models discussed in 


this paper. The model isochrones are curved but can be 
well approximated with a two temperature fit within the 
2800-5000 K region. The slopes presented in Table p 
are intended as a guide, with a measured slope that is 
sensitive to the temperature region defined in the fit. 
For low mass stars (3200-5000 K), the luminosity slopes 
in the Dartmouth, Pisa, BCAH98-1.9, BHAC-2015, and 
Feiden pre-main sequence tracks are consistent with the 
observations. These models produce similar ages for stars 
at 3800 and 4200 K. On the other hand, the isochronal 
slopes predicted from the DM97, Siess, BCAH98, and 
PARSEC models differ from empirical slopes. This dis¬ 
crepancy leads to differences in expected luminosities 
versus spectral type. The DM97, Siess, and BCAH98 
models all have a steeper slope than is observed, which 
leads to a temperature (mass) dependence in age esti¬ 
mates. As an example, if we normalize our empirical 
isochrone and the BCAH98 isochrone to a young 4200 
K star, a co-eval star of 3400 K would appear 0.3 dex 
more luminous (0.4 dex younger in age) than expected 
from the BCAH98 isochrone. The PARSEC isochrones 
conversely have a shallower slope than is observed, which 
leads to 3400 K stars appearing older than 4200 K stars. 

The new BHAC-2015 and Feiden models have 
isochrones that are similar to the empirical isochrones. 
The inflection point in both BHAC-2015 and Feiden oc¬ 
curs at temperatures ~ 50 — 150 K hotter than in the 
empirical isochrones, which leads to mismatches between 
empirical and model isochrones for very low mass stars 
and brown dwarfs. These two models predict isochrones 
that are close enough to empirical isochrones that the 
differences could be explained if effective temperatures 
below 3500 K are systematically underestimated by 50- 
100 K. 

Given these mismatches between model and observed 
isochrones, even for the latest models for cool tempera¬ 
tures, comparisons between clusters should be made at 
a constant temperature. Age estimates for each model 
are evaluated at 3400, 3800, and 4200 K and are plot¬ 
ted in Fig. [12] for three of the benchmark associations 
and for the Upper Sco OB Association. Because of the 
discrepancies in slope described above, for most models 
the 3400 K stars are 0.2-0.5 dex younger than 3800 K 
and 4200 K stars. The PARSEC models have the reverse 
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TABLE 5 

Consistent Sets of Ages (in Myr) of Young Associations' 1 


Stellar 

Association 

d (pc) 

log L4200/±0 

DM97 
3400 K 

Pisa 
4200 K 

BCAH98 
4200 K 

Feiden 
3800 K 

BHAC-2015 
3800 K 

Li Depl. 6 
Age 

Upper Sco 

145 ± 10“ 

-0.25 ±0.06 

1.7 ±0.2 

5.0 ± 1.1 

10.3 ± 1.9 

4.2 ± 1.0 

4.1 ± 1.0 

- 

e Cha 

111 ± 6“ 

-0.28 ±0.06 

1.9 ±0.3 

5.6 ± 1.4 

11.3 ±2.1 

4.7 ± 1.2 

4.7 ±1.2 

- 

77 Cha 

94 ±5“ 

-0.32 ±0.05 

2.1 ±0.3 

6.6 ± 1.3 

12.7 ±2.1 

5.5 ± 1.3 

5.5 ± 1.3 

- 

TWA 

(57) 

-0.40 ±0.04 

2.7 ±0.4 

9.0 ± 1.5 

16.5 ±2.1 

7.9 ± 1.4 

8.1 ± 1.4 

- 

BPMG 

(30) 

-0.65 ±0.04 

6.1 ±0.7 

24 ±4 

36 ± 5 

23 ±4 

22 ±4 

24 ±4 

Tuc Hor 

(51) 

-0.72 ±0.02 

7.7 ±0.5 

32 ±3 

44 ±3 

30 ±3 

29 ±3 

40 ±4 


“For models that produce a younger, intermediate, and old age, plus the recent BHAC-2015 and Feiden models. 


Age assessed at the listed temperature because of temperature dependence in age estimates. 
fc Li depletion boundary ages, see text for details. 

“Estimated uncertainty in median cluster distances. 


dependence, with the 3400 K stars older than the 4200 
K stars. 

The brown dwarf ages are not plotted. For the now- 
outdated BCAH98 models, ages of 3000 K stars are an¬ 
other 0.3 dex younger than the 3400 K stars. On the 
other hand, for PARSEC models ages of 3000 K stars 
are typically 0.1-0.3 dex older than 3400 K stars. How¬ 
ever, the differences between the brown dwarf loci of Up¬ 
per Sco OB Association and the Tuc Hor Moving Group 
demonstrates that the shape and slope of the brown 
dwarf isochrone evolves with time. Fitting a single slope 
to the brown dwarf isochrones without respect to age is 
an overly simplistic approach to characterizing the loca¬ 
tion of brown dwarfs in the HR diagram. The DM97 
models reproduce the evolution in the brown dwarf iso- 
chones especially well, however the absolute ages are 
anomalously young. 

5. DISCUSSION 

Converting the luminosity proxy into a consistent and 
accurate tool for age measurements requires solving two 
problems: ( 1 ) establishing whether the ages inferred from 
luminosities are consistent across a temperature range, 
and ( 2 ) comparing the model-dependent age estimates 
with independent age estimates to establish an absolute 
and correct age scale. 

The empirical isochrone fits presented here lead to a 
self-consistent luminosity ordering for nearby young as¬ 
sociations. Absolute ages are then estimated from com¬ 
parisons to isochrones for model tracks, but these com¬ 
parisons lead to a wide range of age estimates for different 
model tracks and within the same models for stars of dif¬ 
ferent temperatures. If uncorrected, these discrepancies 
lead to significant spectral type (mass) dependences in 
ages of young stars. 

The temperature dependency of the inferred ages 
arises from the mismatch between empirical and the¬ 
oretical slopes in logL vs temperature. The source 
for this mismatch may be the evolutionary models, al¬ 
though spectral type/temperature mismatches between 
different methodologies and uncertainties in the spectral 
type/temperature conversion also likely contribute to 
these discrepancies. The BHAC-2015 models improved 
the prescription for convection, while the Feiden models 
improve upon the thermal structure in the atmosphere. 
Both the BHAC-2015 and Feiden models significantly im¬ 


prove the match between predicted and observed stellar 
loci. Alternately, the PARSEC models apply an ad hoc 
change to the T-r relationship to fit the color constraints 
of old clusters, and similar tweaks may be able to repro¬ 
duce loci of low mass young stars. 

Ultimately, despite recent advances, some model im¬ 
provements are still necessary so that the evolution of 
radius and temperature better reflects the observational 
constraints of young stellar associations. Missing physics 
in the evolutionary models, including magnetic fields and 
spots, may lead to mass dependent discrepancies or abso¬ 
lute errors in age dating young clusters. Age comparisons 
between clusters are self-consistent only when compared 
over specific temperature ranges where the model and 
empirical isochrones match. The BHAC-2015 and Feiden 
models currently offer the best set of pre-nrain sequence 
tracks for relative age estimates, although mismatches at 
< 3500 K still lead to temperature-dependent ages. The 
empirical isochrones derived in this paper offer a method 
to place stars of different temperatures on the same age 
scale, albeit with some uncertainty introduced by the 
possible temporal evolution of the isochronal slopes. 

Table [5] lists sets of ages for each cluster based on a 
consistent application of a single model family at a single 
stellar temperature. The ages are selected from models 
that produce young, intermediate, and old ages relative 
to the average model. Uncertainties in these ages are 
calculated from the statistical uncertainties in the fits, 
as listed in Table , and an estimate for the systematic 
uncertainty in the median distance of the 77 Cha, e Cha, 
and Upper Sco OB Associations. 

The reasonableness of these isochronal ages are 
testable by comparison to Li depletion boundary ages 10 . 
The Tuc-Hor Moving Group has a Li depletion boundary 
age of 41 ± 2 Myr from the BCAH98 evolutio nary mo d- 
els and 38 ± 2 Myr from the DM97 models dKraus et alJ 
I2014T) . so we adopt 40 ± 4 Myr. The j3 Pic Moving 
Grou p has Li depletion bou ndary ages of 21 ± 4 Myr 
from (|Binks fe Jeffriesj |201 4D a nd 26 ± 3 Myr from the 
Feiden models dMalo et al.ll2014all . We adopt an Li de¬ 
pletion boundary age of 24 ± 4 Myr for the j3 Pic Moving 
Group. 

9 The luminosity errors in dex are asymmetric in linear space. 
The listed errors are the median error in Myr. 

10 The Li depletion boundary ages may be systematically over- 
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Fig. 11.— The empirical HR diagrams, along with two-temperature regime linear fits, compared to isochrones from pre-main sequence 
evolution models. The points for rj Cha and e Cha Association members are shifted to be consistent with the isochrone for TW Hya, while 
the (3 Pic Moving Group members are shifted to be consistent with Tuc-Hor Moving Group members. As in Figure 8, the black lines show 
isochrones of log age = 6.0, 6.5, 7.0, 7.5, and 8.0, with 6.0, 7.0, and 8.0 in bold. 
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Fig. 12.— Age estimates of Tuc Hor Moving Group, f$ Pic Moving Group, the TW Hya Association, and Upper Sco OB Association, as 
calculated from eight different models for 4200 K, 3800 K and 3400 K stars (from left to right within each model). The horizontal hatched 
yellow region shows the Li depletion bound ary ages of 24 ± 4 and 40 =b 4 Myr 6 for the (3 Pic and Tuc-Hor Moving Groups (Binks &: Jeffries 
120141 ; IMalo et al.H2014al ; I Kraus et al .1120141) . 


Figure Q2] and Table [5] compares the Li depletion 
boundary ages to the isochronal ages estimated here from 
the different model tracks. No single model of low mass 
stars reproduces the Li depletion boundary age for both 
the Tuc-Hor and (3 Pic Moving Groups. Most models 
underpredict the age of the Tuc-Hor Moving Group by 
~ 10 Myr (or 0.1 dex). On the other hand, most mod¬ 
els reproduce the age of the (3 Pic Moving Group. The 
DM97 models yield especially young ages, likely because 
the isochrones of low-mass stars have temperatures that 
are initially too high (see description in §2.5). Both the 
Feiden and BHAC-2015 models have similar isochrones 
that yield remarkably similar ages. 

When applied consistently, the Li depletion boundary 
ages for the f3 Pic Moving Group implies a ~ 5 Myr old 
age for low mass stars in the Upper Sco OB Association, 
younger than the 11 Myr a ge measured from intermedi¬ 
ate and high mass stars bv IPecaut et al.l (120121) . The 11 
Myr age of Upper Sco would require ages of ~ 40 and 55 
Myr for the (3 Pic and Tuc-Hor Moving Groups. Either 
the intermediate mass or low mass age (based on Li de¬ 
pletion boundary ages of older associations) of the Upper 
Sco OB Association is incorrect, or the identification and 
evaluation of the low mass and/or higher mass stars is 

estimated by a few Myr if the radius dis persion within a cluster 
decre ases the timescales for Li depletion ISorners &; Pinsonneaultl 
120141) . On the other hand, inclusion of magnetic spots may lead to 
30% and 20% underestimates of the ages of the /3 Pic and Tuc-Hor 
Moving Groups, respectively (Jackson & Jeffries 2014). 


biased in some way. 

Errors in ages may be caused by either errors in the 
birthline or in the contraction rates. A bulk subtraction 
of 5 Myr to intermediate mass isochrones would shift the 
intermediate mass isochrones sufficiently for Upper Sco 
and would not cause any significant discrepancy with the 
Li depletion age of the (3 Pic Moving Group. The age 
estimate for the (3 Pic Moving Group from late F and 
G stars dMamaiek fe Bel3 12014T) would also not change 
significantly because these stars are preferentially lower 
mass than the intermediate/high mass stars that led to 
the old age estimate for the Upper Sco OB Association. 
An addition of 5 Myr to low-mass isochrones would also 
lead to consistent ages, but would encounter severe prob¬ 
lems by forcing the youngest known clusters to have ages 
of 5 Myr. Alternately, the contraction rates for either 
low-mass or intermediate mass stars may be in error. For 
low mass stars, the inclusion of magnetic fields and spots 
in evolutionary models will affect contraction rates. For 
intermediate and high mass stars, the inclusion of rota¬ 
tion leads to older ages for intermediate and high-mass 
stars. 

Observational errors relative to the benchmark associ¬ 
ations may also lead to an underestimate of the age. A 
bulk increase in temperatures by 250 K would lead to an 
age of 10 Myr, but such a large systematic difference in 
spectral typing is unlikely. Alternately, the luminosities 
could be overestimated by 0.2 dex if the distances are bi- 
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ased and overestimated by 25%. Assuming Ay = 0 mag, 
which disagrees with the colors of Upper Sco stars, would 
lead to only a 0.07 dex increase in the luminosities. 

6 . CONCLUSIONS 

In this paper, we calculate empirical isochrones in HR 
diagrams for low mass stars and brown dwarfs in young 
(5 — 40 Myr) nearby associations. The sample consists of 
stars that have accurate distances, low extinction, and 
are not affected by the presence of disks or accretion. 
The simple formulation of the empirical isochrone as¬ 
sumes that two slopes, one applied to stars with tempera¬ 
tures 3200/3400-5000 K and one for stars between 2800- 
3200/3400 K, applies to the stellar locus of all young as¬ 
sociations. This choice leads to reasonable fits to the five 
young benchmark clusters and to the stellar and brown 
dwarf locus of the Upper Sco OB Association. 

A successful theory of pre-main sequence evolution 
should match these empirical fits. The Dartmouth, Pisa, 
and BCAH98-1.9 pre-main sequence isochrones between 
3200-5500 are consistent with the measured slopes. The 
slopes from DM97, BCAH98, and Siess isochrones are 
steeper than observed, so that very low mass stars will be 
assessed as younger than solar mass stars. The PARSEC 
isochrones have the opposite behavior, so that solar mass 
stars are assessed as younger than very low-mass stars. 
The PARSEC and BCAH98 models typically lead to the 
oldest ages, while the ages obtained from the DM97 mod¬ 
els are anomalously young. Advances in prescriptions for 
convection and the thermal structure of the atmosphere 
in the BHAC-2015 and Feiden models, respectively, have 
improved the match to empirical stellar loci at > 3500 
K and match the decrease in luminosity at < 3500 K, 
though with a temperature offset relative to established 
spectral type-temperature conversions. 

Applying the empirical isochrone to associations pro¬ 
vides a statistically robust age ordering, which is inde¬ 
pendent of spectral type. Our age ordering is based 
on the 3200-5000 K region. The assumption that the 
2800-3200 K isochrones have the same slope for both 
the young and old associations is likely overly simplis¬ 
tic. Steeper slopes at younger ages, especially for the 
Upper Sco OB Association, instead indicates some evo¬ 
lution in the brown dwarf locus. Because the relationship 
between luminosity and temperature for brown dwarfs is 
steep, any use of brown dwarfs in age estimates requires 
accurate temperatures and should be used with caution. 

Fitting these isochrones to the six associations dis¬ 
cussed in this paper leads to a luminosity/age ordering 
mostly as expected from previous results, with the Up¬ 
per Sco (L 4200/^0 = —0.23, ?] Cha (—0.28), and e Cha 
(—0.32), Associations roughly co-eval, followed by the 
TW Hya Association (—0.40) and the /3 Pic (—0.65 and 
Tuc-Hor (—0.72) Moving Groups. For any given asso¬ 
ciation, age estimates of different models evaluated at 
different temperatures span 0.8 dex. As a consequence, 
comparisons of ages between different associations should 
consistently apply a single luminosity-age scaling. The 
BHAC-2015 and Feiden pre-main sequence models yield 
very similar ages for low-mass stars. Ages can be adopted 
to be consistent with the Li depletion boundary age of 
the /3 Pic Moving Group, but these models imply a ~ 5 
Myr age for the Upper Sco OB assocation for low mass 
stars, younger than the 11 Myr age obtained from inter¬ 


mediate mass and high mass stars stars. 

We thank the referee, Cameron Bell, for a careful read 
of the manuscript, which led to substantial improvements 
in the robustness and clarity of results and in the self- 
consistency within the paper. GJH thanks Adam Kraus 
for discussions on preparing the paper and Jessy Jose for 
valuable comments on the paper. GJH is also grateful for 
the role of the 2014 Oort Workshop on Episodic Accre¬ 
tion, jointly organized by Ewine van Dishoeck and Neal 
Evans, in motivating this analysis, and discussions with 
Lee Hartmann and Isabelle Baraffe at the workshop on 
early ideas for this paper. We appreciate Gregory Feiden 
sharing the new Dartmouth tracks prior to their publi¬ 
cation and for discussion of those tracks and magnetic 
fields. GJH is supported by a Youth Qianren grant and 
general grant # 11473005 awarded by the National Sci¬ 
ence Foundation of China. This research has made use 
of the VizieR catalogue access tool, CDS, Strasbourg, 
France. 

7. APPENDIX A: TARGET PHOTOMETRY AND 
BOLOMETRIC CORRECTIONS 

In §2.2-2.3, we describe the collection and use of pho¬ 
tometry to identify disk sources, to tweak bolometric 
corrections calculated from models, and to measure ex¬ 
tinctions for members of the Upper Sco OB Association. 
This appendix describes how the photometry is imple¬ 
mented in this paper. 

Filter curves and zero points are obtained from the 
Virtual Observatory Filter Profile Service 11 . The BT- 
Settl models were convolved with the filter transmission 
profiles in photon space, which may lead to minor offset s 
relative to integrations in energy space ( Cohen) 12003 ). 
The number of counts is then converted into a magnitude 
from the filter zero point function. Synthetic photometry 
is then calculated from the zero-point flux for the filter 
and the central wavelength of the filter. The listed cen¬ 
tral wavelengths are the average wavelength of the filter, 
calculated by integrating the transmission over the filter 
bandpass in wavelength space. 

Figure fl3l shows temperature-color plots that compare 
the observed photometry to synthetic photometry from 
the B T-Sett l models and , w hen possible, to photome¬ 
try of iPecaut fc Mamaieki (120131) . In cases where the 
observed and synthetic colors (always relative to J) are 
mismatched, the model fluxes at the corresponding wave¬ 
length are corrected to reproduce the observed colors. 
The J-band flux is always kept constant, so all correc¬ 
tions are relative to the J-band. The bolometric correc¬ 
tions are then calculated from the synthetic J-band pho¬ 
tometry and the integrated emission in the full model 
spectrum. 

Finally, the absolute bolometric magnitude Mboi is 
calculated from the distance-corrected J-band magni¬ 
tude, Mj, and the bolometric correction, BC, M^oi = 
Mj + BC. The zero point for Mm (3.05 5 x 10 35 erg 
cm -2 s _1 ) and the assumed Lq = 3.83 x 10 33 erg cm -2 
s _1 are obtained from the 1999 IAU Zero Point Scale, as 
discussed at the webpage of Eric Mamajek listing Basic 
Astronomical Data for the Sun 12 . 

11 http://voservices.net/filter/ 

12 https://sites.google.com/site/mamajeksstarnotes/basic- 
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Fig. 13.— Color-temperature plots for members of young (red circles) and old (blue squares) nearby associations. All colors include J 
to evaluate bolometric correc tions from the J band. Sy nthetic colors from BT-Settl models are plotted as the black line, while colors for 
young (5—30 Myr) stars from [Pecaut &; Mamaiekl <|2013T ) are plotted in green. In cases where the observed colors deviate from the model 
colors, the colors adopted for the bolometric correction calculations are shown as a dashed purple line. Objects with disks, identified from 
a K-W2 or K-W3 excess (Fig.0, are shown here as the small filled circles or squares surrounded by an outer circle or square, 
astronomical-data-for-the-sun 
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